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Abstract 

The alternative to the standard formulation of the quark-parton model 
is proposed. Our relativistically covariant approach is based on the solu- 
tion of the master equations relating the structure and distribution func- 
tions, which consistently takes into account the intrinsic quark motion. It 
is suggested, that the intrinsic quark motion can substantially reduce the 
structure function gi. Simultaneously it is suggested, that the combined 
analysis of the polarized and unpolarized data can give an information 
about the effective masses and intrinsic motion of the quarks inside the 
nucleon. 



1 Introduction 

Measuring of the structure functions is an unique tool for the study of the 
nucleon internal structure - together with the quark-parton model (QPM) giving 
the relations between the structure functions and the parton distributions, which 
represent a final, detailed picture of the nucleon. In this sense, these relations, 
obtained under definite assumptions, are extraordinary important, since the 
distribution functions themselves are not directly measurable. At the same time, 
the standard, simple formulas relating the structure and distribution functions 
are ordinarily considered so self-evident, that in some statements, the both are 
identified. 

The experiments dedicated to the Deep Inelastic Scattering (DIS), are ori- 
ented to the measuring either unpolarized or polarized structure functions. The 
results on the unpolarized functions are well compatible with our expectations 
based on the QPM and QCD, but the situation for the polarized functions is 
much more complicated. Until now, it is not well understood, why the integral 
of the proton spin structure function gi is substantially less, than expected from 
very natural assumption, that the nucleon spin is created by the valence quarks. 
Presently, there is a strong tendency to explain the missing part of the nucleon 
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spin as a contribution of the gluons. It has been also suggested, that the quark 
orbital momentum can play some role as well. Nevertheless, a consistent ex- 
planation of the underlying mechanisms is still missing. During the last years, 
the hundreds of papers have been devoted to the nucleon spin structure, for the 
present status see e.g. 0,0, the comprehensive overviews 0,j| and citations 
therein. 

In the present paper we summarize and update our discussion started in 0- 
, where we have shown, that the standard formulation of the QPM, conceptu- 
ally firmly connected with the infinite momentum frame (IMF), oversimplifies 
the parton kinematics. In we demonstrated that the effect of oversimplified 
kinematics in IMF can have an impact particularly on the spin structure func- 
tion gi, or more exactly, it can substantially modify the relation between the 
distribution and spin structure functions. 

The paper is organized as follows. In the following section the basic kinemat- 
ical quantities related to the DIS are introduced and particularly the meaning 
of variable xb is discussed. In the Sec. |3|we consider proton as an idealized sys- 
tem of the quasifree, massive partons with the four-momenta on the mass shell. 
In the covariant formulation we deduce the relations between the structure and 
distribution functions for the unpolarized and polarized case. At the same time, 
the obtained relations are compared with those derived in the standard, IMF 
approach. In Sec. H., using the results obtained in the previous sections, we 
propose a more realistic model of the proton, in which the internal motion of 
quarks is consistently taken into account. In contradistinction to the standard 
treatment based on the QCD evolution of the distribution functions in depen- 
dence on Q 2 , our model rather aims at describing the part of distribution and 
structure function, which is not calculable in terms of the pertubative QCD. 
In the Sec. the results of the model on the polarized and unpolarized pro- 
ton structure functions are compared with the experimental data and some free 
parameters are fixed. Some additional comments on the model and obtained 
results are done in the Sec. The last section is devoted to the summary and 
concluding remarks. 



2 Kinematics 

First of all let us recall some basic notions used in the description of DIS and the 
interpretation of the experimental data on the basis of the QPM. The process 
is usually described (see Fig. by the variables 

q 2 ^-Q 2 = (k-k>) 2 , xb = $t. (2.1) 

2Pq 

As a rule, lepton mass is neglected, i.e. k 2 — k' 2 — 0. Important assumption 
of the QPM is that the struck parton remains on-shell, that implies 

q 2 + 2pq = 0. (2.2) 
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Figure 1 : Diagram describing DIS as a one photon exchange between the charged 
lepton and parton. 



Bjorken scaling variable xb can be interpreted as the fraction of the nucleon mo- 
mentum carried by the parton in the nucleon infinite momentum frame (IMF). 
The motivation of this statement can be explained as follows. Let us denote 

p(lab) = (po,Pi,P2,P3), P(lab) = (M, 0, 0, 0), q(lab) = (q Q ,q 1 ,q 2 ,q 3 ) 

(2.3) 

the momenta of the parton, nucleon and exchanged photon in the nucleon rest 
frame (LAB). The Lorentz boost to the IMF (in the direction of collision axis) 
gives 

p(inf) = bo,pi,p 2 ,P3), P(inf) = (Pq, P{, 0, 0), g(inf) = (q' , q[, q 2 , q 3 ), 

(2.4) 

where for (3 — ■> — 1 



P0=Pi=7(Po+Pi), P^ = P[ = jM, 7 =l/Vl-/3 2 - (2.5) 
If we denote 

P^~ P{ M ' 1 Dj 

then one can write 

p(inf)=arP(inf) + (0,0,p a ,p3). (2.7) 
Now let the lepton has initial momentum k{lab) = (kg, — feo, 0, 0). If we denote 



v = k — k' and qi = qx, then q^ < and from Eqs. (2.1), fl2.2| ) it follows 

pq p v+ | q L | pi p T q T , . 

XB = Pq- = ~Mv W (2 ' 8) 
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where pr, qr are the parton and photon transversal momenta. Obviously 
k' 2 = (k — qf = k 2 + q 2 - 2k v + 2k Q | q L | = 0, 



q L \ Q 2 M 

= 1 + 7T, = 1 + T XB - 

V ZKqV Kq 



Using this relation the Eq. (^8) can be modified 



Po+Pi . Pi Prqr 

XB = — n — h t Xb ~ ~T7f — 

M k Mv 



(2.9) 



(2.10) 



Prqr 
Xb = x — — cos ip, 



therefore, if the lepton energy is sufficiently high, so pi/ko « 0, one can write 

MV 

where ip is the angle between the parton and photon momenta in the transversal 
plane. 

So, if parton transversal momenta are neglected, xb really represents fraction 
of momentum (2.6) and (2.7). In a higher approximation the experimentally 
measured xb, being an integral over ip, is effectively smeared with respect to 
the fraction x— which is not correlated with ip. An estimation of the second 
term in the last equation can be done as follows. Because 



q 



Q 2 



AM 2 



Q 



~ B l 



then Eqs. (2.9), (2.12) give 



9t / / q \ flL\ z I [AM 2 M 2 \ 2 im ^yi-i'u 



v \ v 



therefore, for M/ko k we obtain 

pq Pa + Pi 2p T x B 



Mv M 



■ cos ip 



(2.12) 



2M 2Mx b 
\Q 2 kl) X *-l^ XB< -[W' (2 ' 13) 



(2.14) 



and 



which gives 



2p T x B 
xb = x ;=- cos <p. 



x -. 2p r p 
— = 1 H == cos w. 



(2.15) 



(2.16) 
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i.e. x = xb for Q 2 — > oo. Therefore can be at sufficiently high Q 2 considered 
as a good approximation of x (and vice versa) . Accordingly, in the next sections 
we shall frequently use the approximation 

pq_ Po+Pl , , 

Mv M ' [ U) 

which will considerably simplify calculation of some integrals. The effect of this 
approximation is well under control, all relevant relations could be calculated 
without it, but at a price, that the sense of some obtained expressions would be 
quite apparent only after a numeric calculation. 

Now, if we assume parton phase space is spherical (in LAB) and an idealized 
scenario in which the parton has a mass m 2 = p 2 , — p\ — p\ — p 2 , then further 
relations can be obtained. 

1) vari able x 



From Eq. (2J5) and the condition x < 1 it can be shown 

,2 



*>£r, (2-18) 



'Pl +P2+P3 SPm 



p% < M 2 (x 



2M ' ri ~ "' v ~ M 2 
Obviously, the highest value of p\ is reached if pt = and 

^/p\ +m 2 +px 

x = V7 = 1 

M 



){l-x). (2.19) 



(2.20) 



which gives 

Pl max — Pm 

Then spherical symmetry implies 



M 2 - m 2 
2M 



pI+P 2 2 +pI <Pm, 



(2.21) 



(2.22) 



i.e. the first relation in (2.19) is proved. Apparently, the minimal value of x 
is reached for pi = —p m an d Pt = 0. After inserting to (2.6) one gets ( 2.18Q . 
Finally, the relation (2^6) implies 



M 2 x 2 -m 2 -Pt 
2Mx 



which, inserted to modified relation (2.22) 



Pi+Pt< 



( M 2 -m 2 s 2 
\ 2M 



(2.23) 



(2.24) 
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after some computation gives the second relation in (2.19). 

2) variable xb 
Let us express is in the LAB 

pq p v-pq If I 2 p q\ 



and estimate its minimal value. With the use of ( priii we obtain 



x B > ^ Wrn 2 +P 2 m -p m Jl + ^-x%) . (2.26) 



Since 



and 



AM 2 9 2M 2 , 

< 1 + -^-4 (2.27) 



^(>^?+^-ft»)=^, (2-28) 



relation (2.26) can be rewritten 



m 2 2Mp m 2 m 2 2Mp m m 4 m 2 ( 2p m m 2 , 
^^^-^^^M2--Q^MJ = M?(, 1 - M-Q^l- (2 - 29) 



i.e. for m 2 Q 2 lower limit of coincides with the limit (2.18) 



3 Idealized scenario: quasifree partons on mass 
shell 

In this section we imagine the partons as a gas (or a mixture of gases) of quasifree 
particles filling up the nucleon volume. The prefix quasi here means that the 
partons bounded inside the nucleon behave at the interaction with the external 
photon probing the nucleon as free particles having the momenta on mass shell. 

3.1 Deconvolution of the distribution function 

Let us suppose F(x) is the distribution fu nct ion of some sort of partons given 



in terms of variable x according to Eq. (2J3) and these partons are assumed 
to have the mass m. If the spherical symmetry is assumed in the nucleon rest 
system and G{po)<fip is the number of partons in the element of the phase space, 
then the distribution function F(x) can be expressed as the convolution 

F{x) = J 5 {^f 1 - xj G( P o)d 3 p, po = ^/m 2 +p 2 +p 2 +p 2 . (3.1) 
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Using the set of integral variables h,po,<p instead of pi 1 p2,p% 



Pi = h, P2 — \Jp\~ m 2 — h 2 sin tp, p 3 = J p 2 , — m 2 — h 2 cos tp, (3.2) 



the integral (3.1) can rewritten 



F(x) = 2n J ™ J + " 5 - x\ G(p ) Po dhdp , H = JpJ 



mr. 
(3.3) 



First, let us calculate inner integral within limits ±H depending on pq. For given 
x and po there contributes only h for which 

p Q + h = Mx, (3.4) 

but simultaneously h must be inside the limits 



p 2 - m 2 < h < Jp 2 - m 2 (3.5) 



which means, that for 



or equivalently for 



Pa + y pi - m 2 < Mx (3.6) 



Mx m 2 



considered integral gives zero. For pq > £, when the both conditions (p. 4), (3.5) 
are compatible for some value h, the integral can be evaluated 

H Jpo + h 



H 



5 ( - x ) G( Po )p dh = MG(p )p . (3.8) 



Therefore the integral (3.3) can be expressed 

F(x) = 2nM { ^ G(p )podp Q . (3.9) 



Let us note, the equation similar to this appears already in || but with the 
structure function F%{x) instead of the distribution one. We shall deal with 
the F2 in the next section, where it will be shown, that the corresponding 
relation is more complicated. For a comparison see also where on the place 
of G(po) the statistical distribution characterized by some temperature and 
chemical potential is used. 

Next, from the relation (3.7) we can express 1 as a function £ 



a±^e^n 2 

x ± = . (3.10) 
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Using the relations ( 2.18Q , ( |3.7| ) one can easily check 



m m 
1 > X-l > — > x- > 



M 



M 2 ' 



M 2 + m 2 

£ max = 2M > e > m. 



First let us insert x+ into (3.9) 



M 



2nM 



G(po)p dp . 



Differentiation in respect to £ gives 



G(0 



1 ^K + VC 2 -m 2 \ / 1 



:F' 



1 



Now we integrate the density G(po) over angular variables obtainin 



P{po)dp = / G(p a )d 3 p = AnG{p )p \J pi - m 2 dp 



and after inserting into (3.13) we get 



P(po)d Po = -2F 



, ( Po + VPu " ™ 2 \ Po + V^o ™ 2 d Po 



71/ 



M 



M 



Second root x_ gives very similar result 

P(po)dp = +2F 
From the definition 



, I Po ~ y/Po m2 1 Pa ~ \Jp\- ™ 2 dp 



M 



M 



M 



x± 



_Po± y/pl - m 2 



M 



the useful relations easily follow 

2 Po 



X+X- 



m 
M 2 " 



X + + X- = 



iVpT 



M 



dpp 
M 



1 m 2 
-M- T7IZT) dx ±, 



dx+ 
x+ 



dx. 

X- 



Now, the equations ( 3.15 ), ( 3.16 ) can be joined 

2 



P(Po) 



T^F'(x ± )x ± . 




Figure 2: Example of the function obeying Eqs. ( 3.23 ), ( 3.2^ ). 



How to understand the two different partial intervals (3.11) of x give indepen- 
dently the complete distribution P(po) in Eq. (3.20)? It is due to the fact that 
e.g. X- represents in the integral (pTTI) the region 



J pi +Pr + m 2 +pi m 
— — = 1 = X- < — . 

M ~ M 



given by the paraboloid 



p\ < 2m \pi 



Pi<0, 



(3.21) 



(3.22) 



containing complete information about G(po) which is spherically symmetric. 
The similar argument is valid for x+ representing the rest of sphere. The Eqs. 
( |3.15 ), (3.16) imply the similarity of F{x) in both intervals 



F'{x+)x+ 
F'{x_)x- 



(3.23) 



which with the use of second relation ( 3.1 9| ) can be easily shown to be equivalent 
to 



F(x + )=F(x_). 



(3.24) 



The relation ( p.20[ ) implies the distribution function F(x) should be increasing 
for (ra/M) 2 < x < m/M and decreasing for m/M < x < 1 e.g. as shown in 
Fig. ||. Now let us calculate the following integrals. 
The total number N of partons: 



N = 



-Em a, 



P(po)dpo = - 



m/M 



F'(x + )(x+ 



M 2 x. 



-)dx-f 



(3.25) 
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= / F' (x+)x + dx+ + / F' (x+)x-dx + . 

Jm/M Jm/M 

The last integral can be modified with the use of ( 3.1 9| ) , ( |3.23 ) 

r 1 r 1 dr^ rm 2 /M 3 

/ F'(x+)x-dx+ = - F'(x-)x 2 _^±= F'{x~)x-dx-. 

Jm/M Jm/M x + Jm/M 

(3.26) 



Then integration by parts gives 



N=- [ F'{x)xdx= [ F{x)dx. (3.27) 

Jm 2 /M 2 Jm 2 /M 2 

The total energy E of partons: 

E = \ P(pa)Podp = F'(x+)(x + - — )~z-(x + + x_)dx + = 

Jm Jm/M M x + 1 



(3.28) 



= ~ / F'{x + ){x\~xl)dx + . 

L Jm/M 



A similar procedure as for N then gives the result 

M I' 1 f 1 
E= / F'(x)x 2 dx = M F(x)xdx. (3.29) 

2 Jm 2 /M 2 Jm 2 /M 2 

Therefore, the both descriptions based cither on the IMF variable x or the 
parton energy po in the LAB give the consistent results on the total number 
of partons and the fraction of energy carried by the partons. Let us remark, a 
model based on the spherically symmetric Gaussian distribution of the parton 
momenta in the hadron rest frame has been recently proposed in fTofl . 

3.2 Structure functions F 2 , Fx 

An important connection between the structure and distribution functions can 
be within QPM derived by a few (equivalent) ways, see e.g. textbooks |rH| - 
[ fL3| . In this section we shall consider the electromagnetic unpolarized structure 
functions assuming quasifrce partons with spin 1/2. The general form of cross 
section for the scattering electron + proton and electron + point like, Dirac 
particle can be written 

da(e- + P ) = -j - K^WcB^M —jy- ^ , (3.30) 

V ; q 4 4^(fcP) 2 - mlKP ' 2fc^) 3 ' V ; 
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da(e + 1) 



K^L^Siip + qf-m')-^ 



<7 4 A\J (kp) 2 — m 2 m 



d 3 k' 



where the electron tensor has the standard form 



K af3 = 2{k a k" 3 + k' a kP + g al3< ^) 



2fc^(27r) 3 ' 

(3.31) 



(3.32) 



and the remaining hadron and lepton tensors W a p, L a p can be written in the 
"reduced" shape 



W, 



a/3 



PgP/3 

M 2 



W 2 - g a pWi, 



(3.33) 



L af 3 = ^p a pi3 - Zgappq- 



(3.34) 



General assumption that the scattering on proton is realized via scattering on 
the partons implies 



da(e~ +p) = J F(£)da(e- +l)d£, 



(3.35) 



where -F(£) is a probabilistic function describing distribution of partons accord- 
ing to some parameter(s) £. Now, if F(£) is substituted by the usual distribution 
function and we assume 



Pa = £,P a ■ 



(3.36) 



then it is obvious, that Eq. ( 3.35 ) after inserting from Eqs. ( |3.30 ) and ( 3.31 ) 
will be satisfied provided that 



W 2 



l 

M 



no 



{2eP a P P - g a p iPq)mP + 1) 2 - ™ 2 R- 

(3.37) 



For simplicity in this equation, and anywhere in this section, the weighting by 
the parton charges is omitted. In fact the Eq. (3.37) is just a master equation 
in fill] (lesson 27, Eq. (27.4)), from which the known relations follow: 



2MW 1 (q 2 ,is) 



F 2 (x) 



xF{x) = F a (a) = vW 2 {q 2 , v), x 



2Mv 
(3.38) 



Here, let us point out, this result is based on the approximation (3.36), which is 
currently accepted in IMF. In fact, relation (3.36) in the covariant formulation 
is equivalent to the assumption, that the partons are static with respect to the 
nucleon, therefore there are suppressed not only the transversal momenta, but 
also longitudinal ones. In the LAB this relation implies £ = m/M, so in the 
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case of our quasifree partons, corresponding distribution function reflects rather 
distribution of the parton masses. 

Before repeating the above procedure for our distribution G(po)d 3 p in LAB, 
one has correctly account for the flux factor corresponding to partons moving 
inside the proton volume with velocity v — p/po- If this velocity has the opposite 
direction to the probing electron, then after passing through the whole subset 
G(po)d 3 p the electron has not still reached backward boundary of the proton, 
where meanwhile the new partons appeared. And on contrary, if the velocity 
of subset has the same direction as the electron, then not all of these partons 
have the same chance to meet this electron. Namely, the partons close to the 
backward boundary are excluded from the game sooner than the electron reaches 
them. Quantitatively, in a subset of partons G(po)d 3 p, the number of partons 
limited by the proton volume and having chance to meet the probing electron 
will be 



dN = (l- v 1 /v e )G( Po )d 3 p, 



(3.39) 



where v\ — px/po is the component of parton velocity in the direction of the 
passing electron, v e ~ — 1 is the electron velocity, if we assume electron momen- 
tum k = (fco, — fco,0,0). Therefore, one can put 



da{e+p) = J (l+ Pl /po)G(po)da(e + q)d 3 p. (3.40) 
We neglect the electron mass, so inserting from Eqs. (3.3C),([T31) gives 



K a(3 W a0 ^K^— [ 
1 2M J 



a gkP [{1+Px/Pq) 



kp 



G( P o)L aP d((p + g) 2 - m z )d 6 p. (3.41) 



The flux factors expressed in the proton rest frame kP = k^M, kp = fco(po +Pi) 
and tensors W a f3,L a p from Eqs. (3. 33), (3. 34) inserted to the last equation give 



= K a f>U G(p )(2 Pa pp - g a0 pq)S((p + q) 2 - ™ 2 )^} 

The last equation can be rewritten 

K a0 (l a(3 - r af} ) = 0, 



(3.42) 



(3.43) 



where l,r are corresponding tensors {...} in the l.h.s. and r.h.s. of Eq. (|3.42|) . 
Since the tensor K a/3 obeys the current conservation 



7 a K af3 = q p K^ = 0, 



(3.44) 
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Eq. (3.43) will be satisfied, if the difference (I — r) a p has the form A{P a qp 
PpQa) + Bq a qp. In this way we get the tensor equation 



W 2 

PaPfi M2 ~~ dapW! + A{P a qi3 + Pf3q a ) + Bq a q p 



(3.45) 



G(p ) 
Po 



(2p a p - g af3 pq)S({p + q) 2 - m 2 )d 3 p, po = \J m 2 +p\ + p% 



for which ( ^3.3€ ) is not required. One can prove, that the tensors W, L in the form 
satisfying ( 3.44 ), would obey Eq. ( 3.45 ) in which A = B = 0. The terms with 
the functions A and B do not contribute to the cross section. Also let us note, 



the veloc ity correction similar to ( 3.3E ) was not used in the Eq. ( 3.37 ) since due 
to (3.36) the all partons in the applied approach have the same velocity as the 
proton i.e. v\ = 0. 

Now the contracting of ( |3.45| ) with tensors g a @ , q a qP , P a P 13 , P a q 13 gives as 
a result the set of four equations 



W 2 ~ 4Wi + 2Mv(A - xB) = 



Mv 



G(p { 



0) r 2 

m — 



2Mxv}5 ( 



pq 

Mv 



x d 3 p, 



(3.46) 



^-W 2 + Wl - 2Mv(A - xB) = ±-J?M [Mxv]8 (£ - x) 



(3.47) 



W 2 - W x + v(2MA + „B) = -L [ \p 2 

Mv J po 



2 1 \Mv 



- x d 3 p, 



(3.48) 



W 2 - Wi + {Mv - 2M 2 x)A - 2MxvB 



(3.49) 



1 fG( Po ) Mxi 

= -TT- / \poMx — 

Mv J po 2 



J \Mv 



x d*p, 



in which the 5— function from the integral ( 3.45| ) is expressed 
S((p + q) 2 -to 2 ) = S(2pq + q 2 ) 



(3.50) 



S 2Mv 



Av ( M. 



\Mv 2Mv 



2Mv \Mv 
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If we define 



(3.51) 



then the solution of the set (3.46)-(3.49) reads 



2MW X 



2Mx + v 



V-i{x) 



M ( m 2 : 



m 2 x 



(3.52) 



T _, .2Mx T , . ,2M 2 x 



ISW2 = X 



2Mx + v 



V-i(x) 



M ( m z 
~ \ M 2 



- 2- 



(3.53) 



... . 6Alx T , . ,6M 2 x 
+Vb(x) + l/i(a;) =- 



2Mi 



V-i(x) 



m 
M 2 



3x' 



(3.54) 



-V (x)2x 1 



Mi 



3Mx\ 



i/ 3 S = 



2M x 



V-i(x) 



m 
M 2 



3a;' 



Mv 



(3.55) 



-y (a:)3a; + Vi(a;) 1 



Mx 



For next discussion we assume v ^> M, then 



vW 2 = F 2 {x) = x 2 F_i(x), MWi = Fi(a;) - (3.56) 

so it is obvious the Callan-Gross relation 2xF\ — F2 holds in thi s approximation. 

In the next step, with the use of the approximation ( 2.17 ), we express the 
integrals fl3.51 ): 



V j (x)= / G(po)(§) S 



j r I Po + Pi 



M 



~ x d 3 p. 



(3.57) 
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This relation with the use of ( 3.1 ) ,( 3.2C ) implies 

M 



where x± is defined in (3.17). The relations (3.58) and (3.18) give 



YM. 



j-k 



( f± 









~2\ J'- fe 






(I- 




2 






2x / 



a; 



(3.58) 



™ ,3,9, 



In the previous section we have shown such functions as Eq. (3.57) obey the 
relation (3.24), which means in particular, that the functions have a maximum 
at xq and vanish for x < Xq. Therefore the same statement is valid also for 
functions F 2 /x 2 and Fi/x from Eq. (|3.56 ) 



F 2 (x+) F 2 {x-) 



Fi{x ± ) = fi(x-) 

x + X- 



(3.60) 



This means that the structure functions of our idealized nucleon have the max- 
imum at xq or higher. 

Further, our considerations have started to move in previous se ction from 
the distribution function F(x) for which we hav e obt ained relation ( 3.2C ). The 
combination of this equation with ( 3.56| ), ( 3.58 ) and ( 3.1 8| ) gives 



1 fF 2 (x)V 



(x 2 +xl), X 



Po 



M 



(3.61) 



(3.62) 



How do we compare the last equation with the standard relation (3.38) for F and 
F2? As we have already told, the standard approach ( 3.37 ) would be exact in 
the case when th e par tons are static with respect to the nucleon, i.e. when x = 
m/M. The Eq. ( 3.45 ) itself is more exact, but in further procedure we assume 
the masses of all the partons in the considered subset being equal. Therefore 
for a comparison let us consider first the extreme scenario when the parton 
distribution functions F(x) and P(po) are [(see Eq. (3.20)] rather narrowly 
peaked around the points xq — m/M and po — m. Then for x ~ xq Eq. ( 3.62 ) 
gives 



lFMx), . FUx) 

- — 2~ (xo + x ) = (3.63) 

2 I X 



from which the second relation (3.36) follows as a limiting case of (3.62) 



x F(x ) w F 2 (x ). 



(3.64) 
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Now, in the case when the distribution functions are broad, the exact validity of 



Eq. (3.37) again requires static partons, therefore the corresponding distribution 
function represents also a spectrum of masses. But then obviously the above 
procedure for a single to can be repeated with spectrum of masses F(xq) giving 



in the result instead of Eq. (3.64) the relation 



x F(x )S(x - x )dx = / F 2 (x )S(x - x )dx , 



which implies 



xF(x) =F 2 (x). 



(3.65) 



(3.66) 



In this sense the standard approach based on Eq. (3.37) can be understood as 
a limiting case (static partons) of that based on Eq. ( [3.45 ) . 



3.3 Spin structure functions g\,g2 



In the previous section t he m aster equa tion ( 3.45 ) has been based on the stan- 
dard symmetric tensors ( |3.3S ) and (3.34) corresponding to the unpolarized DIS. 
After introduction the spin terms into both the tensors (see e.g. 11 , Eqs. (33.9), 
(33.10)) the master equation reads 



W 2 > 



S a MG 1 + {PqS a - SqP° 



M 



(3.67) 



+A{P a q,3 + Pf3q a ) + Bq a q/3 = J G{p a )(2p a p /3 - g a ppq)8({p + q) 2 - to 2 ) 



fp 
Po 



+ie a p\ a q X I H(po)mw' y 8({p + q) 2 - to 2 ) — , 

Pa 



where G and H are related to the polarized quark distributions 
G boHE e ^ T (Po) + ^(Po)), 



(3.68) 



H(po)=J2 e2 A h }(P°)~ h iM) 



(3.69) 



and the spin polarization vectors satisfy 

-1, = = 0. 



(3.70) 
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TheEq. d3~67] ) requires for the spin terms 



S'MGi + (PqS* - SqP*)^ = 

M 2M v 



H(p ) 
Pa 



pq 



w^S[^--x)d"p, (3.71) 



where we use for the S— function the relation ( 3.50 ). 

Now, let us consider first a extremely simple scenario (in LAB) assuming 
the following. 

1) To the function H in Eq. (3.69) only the valence quark term contributes. 

2) Momenta distributions have the same (spherically symmetric) shape for u 
and d quarks 



1 



hd{Po) = ^K{po) = Hpo) 

and both the quarks have the same mass m. 

3) All the three quarks contribute to the proton spin equally 



(3.72) 



1 



1 



K-K= -(hi - hi) EE Ah( Po ) = -h(po) 



p = \lm 2 +pl+pl +P§, 



(3.73) 



i.e. in a first step we ignore constraints due to axial vector current operators 
on the spin contribution from different flavors. Since all the three quarks are 
assumed to give the proton spin, the last equation implies 



3 J Ah(p )d 3 p = 1. 



The combination with ( 3.69] ) gives 



H{ Po ) = 2^Ah( Po ) + -Ah(p ) = Ah{ Po ) 



and 



J H(p )d 3 p=± 



(3.74) 



(3.75) 



(3.76) 



Now, let us assume the proton is polarized in the direction of the collision 
axis (coordinate one), then Eqs. ( |3.70| ),( ^71 ) require for the proton at rest 



S = (0,1,0,0) 
and for the quark with four-momentum p, 

' Pi Po 



Vpo - pi Vpo - pi 



,0,0 



(3.77) 



(3.78) 



More rigorous derivation of this form of the quark polarization vector, which 
is based on the requirement of the relativistic covariance, is done in the next 
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section. The contracting of Eq. (3.71) with P a and S a (or equivalently, simply 
taking a = 0, 1) gives the equations 



qi G 2 



rr> 



H(po) 



Pi 



2Mv J po ^pj~pj 



5 (£h~ x ) d3p ' (3J9) 



(3.80) 



2Mv J p Q ^pl~pl 
In the next step we apply the approximations from the Eqs. (|2.9|) and ( 2.1 7| ) 

PI Po + Pi 



q 1 ~ -v, 



Mv 



M 



(3.81) 



Let us note, the negative sign in the first relation is connected with the choice 
of the lepton beam direction giving the Eq. (2.17). The opposite choice should 
give 



qi +v, 



pq _ Po-Pi 
Mv ~ M 



(3.82) 



and one can check the both alternatives result in the equal pairs G±, G 2l which 
read 

o ( n f H (Po) P0+P1 x fPo+Pi \ ,3 ,„ 
2g 1 (x) = 2M vGi = m / — l = = S [ — — x d*p, (3.83) 



Po y/p 2 - p\ 



A I 



2g 2 (x) ee 2M^G 2 = -m / ( ™^ x ) <Pp. (3.84) 



PO y/p 2 Q - p{ 



Po +Pl 



M 



Let us remark, the integration of Eqs. (3.79) and (3.84) over x gives on r.h.s. 
the integral 



H(p ) Pi 



P° Vp 2 o - Pi 



d z p = 0, 



(3.85) 



which is zero due to spherical symmetry. Therefore in this approach the first 
moment of g 2 (x) is zero as well. Now we shall pay attention particularly to the 
function g\, which can be rewritten 



= if / h (Po) — 
3 J Po 



M po+Pi,(Po+Pi 



5[^ 1 f^-x\d s p, x = ^-. (3.86) 
pa - pi \ M J M 



What our assumptions l)-3) do mean in the language of the standard IMF 
approach? In the previous section we have suggested that our approach is 
equivalent to the standard one [based on the approximation ( 3.36j )], for the static 
quarks described by the distribution function h(po) sharply peaked around m. 
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In such a case the last equation for pg sa m, p\ ~ after combining with (3.69) 
and (l3.ll) gives 



2Si(*) = | 5>2(ftT(po)-*)(po) 



1^ .')</> C-S,S7, 



where /j(x) are corresponding distribution functions in the IMF, so in this 
limiting case our spin equation ( i.8(\ ) is again identical to the standard one, see 
Eq. (33.14) in pi|. The last equation can be in our simplified scenario rewritten 



2si (x) 



h(p )5 



M 



-x)d A p= -f{x) 



F 2 val (x) 

3x 



(3.88) 



This relation says, what our assumptions l)-3) mean in the terms of the IMF 
approach, in particular we obtain 



V i ,,-/ = / gi{x)dx = i J f{x)dx = ^. 



Further, in accordance with (3.57) let us denote 



x dp, 



then ( |3~56| ) and < ^8§j ) give 

, . xV_i(x) 

251(1) = — 



Fimf = g / xF_i(x)dx. 



(3.89) 



(3.90) 



(3.91) 



So, our Eq. (3.86) in the limit case of the static quarks coincides with 
the standard IMF approach, but what this equation implies for the nonstatic 
quarks? Let us calculate the first moment of our gi : 



Xq 



M I p a + pi ( po+pi 



r,„, f j J h{po y^i2L^ 5 [!^-x)d*pdx. 



(3.92) 



Due to the S— function, the square root term in the integral can be rewritten 

(3.93) 



P0+P1 



Mx 



Po - Pi 



2p - Mx 



IMx 
2p~o~ 



Mx\ 
2p~o~) 



-1/2 
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and using Eq. (3.90) the integral is, correspondingly 



lab 



(3.94) 



The integration by parts combined with the relations ( 3.59 ) gives 

J V-js/zix) dx 



2{x/2) 



3+3/2 



3-3/2 W j+3/2 



dx = J Vfa) ( 



- + — 

2 2x 



-J-3/2 



2(x/2) 



3+3/2 



j + 3/2 



-da; 



.7+3/2 Vl + ^o/^ 2 



V {x)2 



3 + 1/2 



3+3/2 



1 + xl/x 2 



2xl/x 3 



{l + xi/xiy 



rdx. 



If we denote t = x\jx 2 and z = 1/(1 + t 2 ) then Eq. ( 3.94 ) can be rewritten as 
r ia6 = ~/ V (x)4t 3 z 5 / 2 f^ ^y^y^dx = I J V Q {x)U s z 2 ^^-dx, 



which implies 



r 



lab 



1 f 1 Ax 2 /x 2 



6U (l+x 2 /x 2 ) 2 



Vo(x)dx. 



Simultaneously, since 



V (x)dx = - [ xVo(x)dx = - [ xVLAx) (- + — dx 
I hi hi V 2 2x J 



V!_ x {x) — + )dx = / V-i(x)xdx, 



the integral (3.91) can be rewritten as 



(3.95) 



6 



Vo(x)dx. 



(3.96) 
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Let us express the last integral as 

Vo(x)dx = I Vo(x)dx + I Vo(x)dx 



and modify the first integral on r.h.s. using substitution y — Xq/x 



rVo(x)dx= fvJ X ^\ X \dy. 

Jx% Jx \y J y 



Now let us recall the general shape of the functions (3.90) obeying Eq. (3.24), 
which implies 

so instead of Eq. ( |3.96 ) one can write 



IMF 

6 



1 fl Vo(x)(^±^) l h, { 3.97 i 



x 



Similar modification of Eq. ( 3.95 ) gives 



lab 

6 



1 t\r,-J H 



M*) o , 2 dx - ( 3 - 98 ) 



Obviously both the integrals are equal for Vq sharply peaked around x = xq, 
but generally, for nonstatic quarks 

Tiab < IW. (3.99) 

Therefore, starting from the one test quark distribution function Vo we get the 
two different results on the first moment of the function g\ - depending on the 
used relation connecting distribution and structure functions. 

What can our result (H||) mean quantitatively? Obviously it will depend 
on the function Vo, on its width. To get some feeling, we use for the Vq the 
following parameterization. According to the Eq. (3.20) for x > xq one can 
write 

xVfa) = -yP(po), PO = Y ( X + \j) ' ( 3 - 10 °) 

Now, for po close to m let us parameterize the energy distribution by 

P(p ) w exp (-Q— ) (3.101) 

satisfying the normalization 

POO 

/ P( Po )dpo = 1. (3.102) 

J m 
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Obviously, the distribution (3.101) means the average quark kinetic energy 
equals to m/a. Inserting (3.101) into ( 3.100| ) gives 



Vo(x) 

Let us note, for \y\ <C 1 



exp 



H 


X 




) 




x 


X 





(1 + y) a « exp(ay), 



therefore if we substitute the exponential function in (3.103) by 

a/2xQ 



exp 



H 


X 




y 




X 


X 





(l-x) 1 



f(x,x ) 



(3.103) 



(3.104) 



Xq < X < 1, 



then the resulting Vb(x) will coincide with ( [3.102 ) in a vicinity of x$, i.e. for 
small kinetic energies, but moreover will obey the global kinematical constraint 
fl3.24| ) outlined in Fig. §. T he rat io of i ntegra ls ( $.98\) and ( |3.97| ) calculated by 
parts with the use of Eqs. ( |3. 103 ) and (3.104) gives 



R s (a,x ) 



lab 



r 



IMF 



4 Ixo x °l x (arctan[a;/a;o] - 7r/4) f(x,x )dx 

Sxo ^ ~ X o/ x2 ) f( X > X 0) dx 



(3.105) 



The results of the numerical computing are plotted in the Fig. ||. What do 
these curves mean? It is obvious, that for static quarks, for which a — > oo and 
R s — > 1 both the approaches are equivalent, as we have already shown. On 
the other hand, it is also apparent, that for nonstatic quarks, with small ratio 
olk, raj (Ekin), both the approaches can differ substantially. 



3.4 Covariant formulation 



Master equation ( 3.67 ) is assembled for quark momenta distributions G, H in the 
nucleon rest frame, but despite of that, the equation is relativist ically covariant. 
Its manifestly covariant form follows immediately from Eq. ( 3.67 ) after the 
substitution 



For moving nucleon we have 

P = (P , i\,0,0), 

and 

pP _ PqPq -PiPi 
M ~ M 



Pi 
Po" 



Pa 
M 



l(P0 - PPl) = PO.iab, 



(3.106) 



(3.107) 



(3.108) 
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a 



Figure 3: Ratio R s plotted for values xo=0.2, 0.02, 0.0002 - in order from top 
to bottom. 



which means, that the phase space of the subset of quarks with po,iab fixed, 
represented by the sphere 



2,2,2 

Pi + Pi + P 3 



P0,lab 



(3.109) 



in the nucleon rest frame, is in a boosted system correc tly rep resented by the 
ellipsoid with the shape defi ned b y the Lorentz transform ( 3.108|) . Let us remark, 
in the same way as the Eq. ( |3.67| ) a similar equation can be obtained and solved 
also for the set of the neutrino structure functions, nevertheless in this paper 
we consider only the electromagnetic ones. 



Further, one can obtain also covariant solution of Eq. ( 3.67 ) for the spin 
functions Gi,G2, but first it is necessary to define correct form of the quark 
polarization vector w. Generally, this vector should be constructed from the 
proton momentum P, proton polarization vector S and the quark momentum 
p: 



APu + BSu + Cp^, 



(3.110) 



where A, B, C are invariant functions of P, S,p. Then contracting of Eq. (3.71) 
with P& , q a and S a gives the equations 



S « MG ' = TEkJ H {w) <- 4M " + B5 « + c ™> 1 (m ~ x ) $■ (3 - 112 » 
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-MGi - vG 2 



2M 

Elimination of G\ , G2 gives 



™ [h(^)(-B + CpS)s(^--x]^. (3.113) 
2Mv J \M J y 1 ' \Mv ) po K ' 



vP/M 
Sq 







(3.114) 



and since P, g, 5* are independent, C must be zero. The remaining invariants 

(3.115) 



A,B follow from Eq. (3.7C), which implies 

A 2 M 2 - B 2 = — 1, ApP + PpS = 

and solution of these equations reads 

pS pP 



A = T - 



B = ±- 



^/( P P) 2 - ( P S) 2 M 2 ' V(P p ) 2 - {pS) 2 M 2 
So the quark polarization vector has the form 

(pP)S^ - (pS)P^ 



±- 



\/(pP) 2 - ( P S) 2 AP 



(3.116) 



(3.117) 



Contributions of both possible solutions [sign +(— ) means that quark spin is 
parallel (antiparallel) to the proton spin in its rest frame] are in our calculation 
taken into account by the difference in Eq. ( 3.69| ). Apparently, for the proton 
r est f rame and polarization S = (0,1,0,0) the last equation is identical to Eq. 
( [3.78 ) . No w, the o btained invariants A, B, C give the spin structure functions 
from Eqs. (|3.11l|) , ( |3.112[) in covariant form 



Gi = 



2M(Pq)(Sq) 



H [Z- 



pP\ ( P P)(Sq) - ( P S)(Pq) 



M J y/(pP)* - ( P S) 2 M 2 



d 3 p 

Po ' 
(3.118) 



Go = 



mM 



2(Pq)(Sq) 



HI** 

M 



P S 



y/(pP) 2 - { P S) 2 AP \Pq 



^. (3.119) 

Po 



Apparently, according to these relations the structure functions can depend also 
on mutual orientation of S and q. Of course, this dependence is more compli- 
cated, apart the factor Sq ahead of the integrals, integration involves also the 
terms pS and pq. This question is being studied and will be discussed in a sep- 
arate paper. Our further considerations will be based on the results obtained 
in the previous section, which follow from Eqs. fl3.118p , ( |3.119| ) applied in the 
proton rest frame for the longitudinal polarization S = (0,1,0,0). Obviously, 
for this case the last two equations are equivalent to Eqs. ( 3.79| ), ( [3.80 ). 

The scheme based on the Eqs. (3.67) and (3.117) with all their implications 
suggested in the previous sections can be a priori valid for quasifree quarks (on 
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mass shell) filling up the nucleon volume. In this sense the scheme represents 
a covariant formulation of the naive QPM. We have shown that Eq. (3.67) 
in which the quark internal motion is consistently taken into account implies 
the relations between the structure and distribution functions different from 
those obtained in the standard procedure relying on the preferred reference 
frame, IMF, which is based on the approximation p^ = xP^. In the covariant 
formulation this approximation is equivalent to the assumption, that the partons 
are static with respect to the nucleon. Of course, this consequence is somewhat 
obscured just in the IMF, where all the relative motion is frozen, since all 
the processes run infinitely slowly - including the passing of the probing lepton 
through the nucleon. Let us remark, the standard relations (e.g. F2 — x e iQi) 
obtained in the naive QPM with static quarks are currently applied even in the 
standard approach based on QCD improved QPM, which is not a consistent 
procedure, since it means that correct dynamics is combined with incorrect 
kinematics. 

In this way we have shown, that the relations between the structure and 
distribution functions can be, at least on the level of the naive QPM, strongly 
modified (particularly for the polarized case) by the parton internal motion. 
This result can be instructive by itself. Let us remark, the impact of the quark 
intrinsic motion on the function gi (x) has been discussed also in some other ap- 
proaches and necessity of the covariant formulation for the spin struc- 
ture functions has been pointed out in |20|]. 



4 Model 

To illustrate, that the scheme suggested in the previous section can be really 
valid for quasifree fermions, let us look at the Fig. ||, where the "structure 
function" of the deuteron measured in quasielastic e~d scattering j2lj] is shown, 
clearly proving the presence of two nucleons in the nucleus. The similarity with 
the Fig. H is apparent. 

Of course, in the case of partons inside the nucleon the situation is much 
more delicate. The interaction among the quarks and gluons is very strong, 
partons themselves are mostly in some shortly living virtual states, is it possible 
to speak about their mass at all? Strictly speaking probably not. The mass 
in the exact sense is well defined only for free particles, whereas the partons 
are never free by definition. However let us try to assume the following. The 
relations obtained within the scheme suggested in the previous sections can be 
used as a good approximation even for the interacting quarks, but provided that 
the term mass of quasifree parton is substituted by the term parton effective 
mass. By this term we mean the mass, which a free parton would have to have to 
interact with the probing photon equally as the real, bounded one. Intuitively, 
this mass should correlate to Q 2 : a lower Q 2 roughly means, that the photon 
" sees" the quark surrounded by some cloud of gluons and quark-antiquark pairs 
as a one particle - by which this photon is absorbed. And on contrary, the 
higher Q 2 should mediate interaction with more "isolated" quark. Moreover, 
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Figure 4: The structure function for quasielastic e d scattering, see text 



we accept that the value of the effective mass can even for a fixed Q 2 fluctuate 
- e.g. in a dependence on the actual QCD process accompanying the photon 
momentum transfer. This means, that the terms in the relations involving the 
mass of quasifree parton xq — m/M will be substituted by their convolution 
with some 'mass distribution' /i : 



f(x,x )—> / p,(xo,Q 2 )f(x,x )dx =F(x,Q 2 ). 



(4.1) 



In the following we shall propose a simple, but sufficiently general model for 
the unknown distributions /i, G, H, in which all the dynamics of the system is 
absorbed. Then, these distributions will be used as an input for the calculating 
of the corresponding structure functions. Construction of the model is based on 
the following assumptions and considerations: 

1) Parton distribution P(e)de representing the number of quarks in the en- 
ergy interval (e, e + de) can be formally expressed: 



(4.2) 



where Tj is a probability that the nucleon is in the state with j partons (quarks 
+ antiquarks) of various flavors, and pj is the corresponding average one-parton 
distribution, which satisfies 



Pj{e)de = 1. 



(4.3) 



2) Nucleon consists of the three quarks and partons (gluons + qq pairs) 
mediating the interaction between them, as sketched in the Fig. |^a, where 
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Figure 5: Nucleon consisting of the valence and sea quarks - with different 
resolutions, see text. 



the individual pictures represent some terms in the sum ( |4.2| ). The flavors and 
spins of all the quarks in each the picture are mutually cancelled, up to the three 
quarks giving additively the corresponding nucleon quantum numbers. These 
three quarks are in the figure marked by black and in our approach are identified 
with the valence quarks. The reason, that such identification is quite sensible, is 
the following. Apparently, the sum (^^) can be split into quark and antiquark 
parts P q (e), Pq(e), then our valence term reads 

Pval(e) = P g (e) - Pg(e), (4.4) 

which in the x— representation exactly corresponds to the current definition of 
the valence quarks. Correspondingly, the unmarked quarks are identified with 
the sea quarks. But both the kinds of quarks have the same energy distributions 



Pj(e) entering the Eq. (4.2), in this sense they are completely e quiv alent. On the 



other hand, it is obvious, that for the valence quarks, in Eq. (4.2) only "black" 
quarks from the figure contribute, therefore if pj is in the first approximation 
assumed to be independent on the flavor, then 

oo 

P val (e) = 3j2r j pM (4.5) 

3=3 

3) The quarks carry only part of the nucleon energy (mass), 

J P{e)ede = c q M, (4.6) 
the rest is carried by the gluons. In the first approximation we shall assume this 



factor is valid also for any term in the sum (4.2), 



j J Pj (e)ede = c q M, (4.7) 
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which in other words means the ratio of the total energies of quarks and glu- 
ons, together constituting the nucleon mass, is the same for all possible states 
sketched in Fig. ||a. 

4) We assume all the quarks in the nucleon state j have approximately the 
same effective mass xq — rrij/M. One can expect, for higher j the parameter 



Xq will drop and so the sum (4.5) can be substituted by the integral 



Pvai(t) = 3 / p v (x Q )p(e,x )dx , 



pv(x )dx a = 1. 



(4.8) 



Obviously, Eq. (4.2) can be with the use of Eq. (4.7) rewritten in a similar way: 



P(e) 



p(x )p(e 1 x )dx , 



where 



f \ ( \ c i M 



e(x ) = / p(e,x )ede. 



(4.9) 



(4.10) 



The physical meaning of the distributions pv,p is the following. The distribu- 
tion p(xq) represents a probability, that the effective mass of the quark, which 
the probing lepton interacts with, is Xq or alternatively, p(xo)dxo is the number 
of quarks in the interval (xq, xq + dxo), which the lepton has chance to interact 
with. On the other hand, the (normalized) distribution py{xo) can be inter- 
preted as a probability, that the exchanging photon " distinguishes" the quarks 
with the effective mass xq - as expressed by the pictures with different granu- 
larity in Fig. |H|a. In this sense, each picture in the Fig. ||a can be labeled by 



some xq, equally as the corresponding term p(e,xo) in the integral (4.8). Obvi- 



ously, at the same time the pv(%o) represents also the distribution of effective 
masses corresponding to the valence quark term. Intuitively, the probability of 
different contributions in Fig. |B|a should depend also on Q 2 (higher Q 2 ^'better 
resolution'), so we expect 



p(x ) -> p(x , Q 2 ), pv(x ) -> pv{xo, Q 2 )- 



(4.11) 



In the next, we shall identify these distributions with that introduced in Eq. 

(0)- 

5) The relations ( 3.56| ), ( 3.59 ) and ( 3.2C ) give the recipe how to obtain the 
structure function F% from a given energy distribution of the partons with some 
fixed value xq and charge e q : 



F 2 (X, X ) = e 2 lfi(x, X ), tp(x, X ) 



1 2tfl(Qe 
x £ 2 + x 



2Vim = ±Mp{e,x ), 



(4.12) 



2N 



where the sign +(— ) in the second relation refers to the region £ < xq (£ > Xq). 
For the application of this procedure to Eqs. (4.8), (4.9) one has to weight the 
contributions integrated over xq by the corresponding (mean) charge squared. 
Apparently, the charge weight of the valence quarks is the constant 



W va l = ~ 



(4.13) 



for the proton and similarly for the neutron, w va i = 2/9. For the sea we assume 
in the first approximation the " equilibrated mixture" of the quarks u : d : s — 
1 : 1 : 1, so 



1 



Then for the nucleon with j quarks we get 
3w vai + (j - 3)w s 

ea 



3(w va l ~ W sea ) 



(4.14) 



(4.15) 



or in terms of xq 



W(X ) = W sea + 3(lV va l - W sea ) 



e(zo) 
c q M ' 



(4.16) 



Therefore, the energy distributions (4.8), (4.9) generate the corresponding struc- 
ture functions: 



F 2v ai{x 1 Q 2 ) = 3w va i / ^ v (x ,Q 2 )'p(x,x )dx , 



(4.17) 



F 2 (x, Q 2 ) = J ^^w sen +3(uw - w sea )j fiv{xQ,Q 2 )(p(x,x )dx . (4.18) 

6) Now, let us pay attention to the spin structure functions. According to the 
concept suggested in the item 2), only valence quarks contribute to the nucleon 
spin. First, we shall consider the spin functions generated by the valence quarks 
with some fixed effective mass xq, then we shall easily proceed to the case with 
the distribution /xy(xo). 



In Sec. 3.3 we have simply assumed all the three valence quarks contribute 



to the proton spin equally [Eq. (3/73)]. On the other hand it is obvious the 



quark symmetry group can impose an extra constraint on the contributions of 
different quark flavors as it follows e.g. from the philosophy of the well known 
Bjorken |2^] and Ellis- Jaffe (2^] sum rules based on the symmetries U(6) and 
SU(3). If we do not assume any particular group of symmetry, then the different 
spin contributions of u- and rf-quarks can be expressed by the free parameter 
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a, < a < 1, having in the notation of Eq. (3.73), e.g. for the proton, the 
following sense 



Ah u ( Po ) = 2ah(p ), Ah d ( Po ) = (1 - 2a)h{p ), 
where h is the valence distribution 

u(p ) = d(p ) = h(p ), / h(p )d 3 p = 1, 



(4.19) 



(4.20) 



which is not, due to different normalization, identical with the distribution p(e), 
but the both are simply related 



p(e) = Aftey/ e 2 — m 2 h(e), 



(4.21) 



in the same way, as the distributions P, G in Eq. ( 3.14 ). 
In the case of proton, there are the particular cases: 

a) a = corresponds to the mutual spin orientation of the three valence quarks 
(s u ,s u ,s d ) = (-1, +1, +1). 

b) a = 1/3 corresponds to the oversimplified scenario studied in Sec. 3.3, assum- 
ing the equal contribution of all the three quarks; (s u , s u , s c i) = (+1/3, +1/3, +1/3). 

c) a = 2/3 corresponds to the non-relativistic SU(6) approach. From the wave 
function 

| p, T) = ^= (±= \{ud + du)u - 2uud) «l -L |(t| + |t) T -2 TTI> 



+_L |(„d-d u ) u )®-L|(u-,LT) T> 



(4.22) 



one can easily show the mean value of the spin carried by the d(u)— quarks is 
— 1/3 (4/3 ), i.e. (s u , s u , Sd) = (+2/3, +2/3, —1/3), which agrees with a — 2/3 in 
Eq. (p~19l) . 

d) a = 1 corresponds to the mutual orientation of the three quarks (s u , s u , Sd) = 
(+1,+1,-1). 

So, the proton spin function H(p ) en tering the equations ( 3.83| ), ( |3.84 ) and 
expressed in terms of the functions ( 4.19| ) reads 

HP(p ) = ^au(p Q ) + 1(1 - 2a)d( Pa ) = 1(1 + 6a)h( Po ). (4.23) 

Assuming the neutron is isospin symmetric, its corresponding spin function will 
be 

H n ( P o) = ~(1 - 2a)u( Po ) + ^ad(p ) = i(4 - 6a)%o), (4.24) 

therefore the corresponding equations for the nucleon spin structure functions 
read 



g?(x,x )=v$ pin ip j (x,x a ), gj(x,x ) = w™ pm ipj(x,x ), j = l,2, 
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«'' 



(l + 6o) 



'spin 



(4 -6a), 



(4.25) 



where, in accordance with Eqs. ( 3.85 ), ( 3.84 ) 



m f h(p ) pa + p! ( pp+pi > 3 
^i{x,x ) = — W 5 — — x)dp, 4.26) 



M*,**) = -~ fhiPo) 



pi 



Pa+Pi 
M 



x d 3 p. 



The function ipi(x,Xo) can be with the use of Eqs. djgg ), ( |3~93| ) expanded 



(4.27) 



2 



V>i(i,afo) = 
Since Eq. ( |3.59 ) implies 

»W> = -/><«) (^) 



E 7H) 3 Va/ 2 W@ • (4-28) 



J+3/2 



one can easily show the sum in Eq. (4.28) gives 
ipi(x,x ) = -x 



xi 



1 v^m I 



Similar manipulation with the function ip 2 gives the result 



ip 2 (x,xo) 



■TO 



1 TO^ e 2 + ^-2.< 



(4.29) 



(4.30) 



(4.31) 



Obviously, for the case with the distribution fxy, the corresponding spin struc- 
ture functions read 



gj{x,Q 2 



w 



sprit 



Hv(x ,Q 2 )ipj(x,xo)dx Q , 3 = 1,2. 



(4.32) 



Let us note, the structure functions F2, Fival, 9i>92 are n °t independent, all of 
them are in the corresponding way generated by the distributions [iy and Vq 
(or, equivalently by p). 

7) Now, to make the construction suggested above applicable for a quanti- 
tative comparison with the experimental data, we have to use some reasonable, 
simple and sufficiently flexible parameterization for the unknown functions \iy 
and Vo- We suggest the following, 
a) Normalized distribution fj,y is assumed in the form 



Vv{x ,Q 2 ) 



r ( r + S + 2 ) rf, \s 

r(r + i)r(« + i) ' x ° {1 ~ Xo) ' 0<x ^ <1 ' 



(4.33) 
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where the Q 2 — dependence is involved in the parameters r, s. 

b) For the function Vq(x)x we shall use the parameterization suggested in Eqs. 

(poiD-(po^) 



Vq{x)x = TC-norm ' f(x, X ), f(x, X ) = 



(i-x) i-±a 



o./2xq 



(4.34) 



where the upper (lower) sign in the first relation refers to the region x > 
xq (x < xq) and 



f(x,x ) 1 



dx 



(4.35) 



which follows e.g. from Eqs. fl3.102| ), ( |3~25| ). Now, apparently one has to 
accept the parameter a ~ m/ (Ekin) depends on Xq as well. Let us consider 
the following. Sequence of the pictures in Fig. |B|a can be understood as the 
pictures of the one and the same nucleon, but "seen with different resolutions" 
as outlined in Fig. Then, it is natural to assume the momentum P of the 
parton from some picture can be obtained from the momenta p\ of n partons in 
a picture more rightward, representing the parton "seen with better resolution" : 



P 



x 



PX- 



Obviously, the mean values satisfy 



(Po) = n (p 



'AO/ 



P) = C 



■corr ■ n (\p\\) , 



Ccorr — 1 i 



(4.36) 



(4.37) 



where c corr equals 0(1) for the extreme case, when the motion of the partons 
in the corresponding subset is completely uncorrelated (correlated). The last 
relations imply the effective masses and kinetic energies obey 



m(P) > n ■ m(p), (E km (P)) < n ■ (E kin (p)) , 



(4.38) 



which means the quantity a is a non-decreasing function of xq. In this moment 
we know nothing more about this function, in the next section we shall show, 
that a reasonable agreement with the experimental data can be obtained with 
the parameterization: 



a(x ) = ai (-ln(xo))" 



(4.39) 



Since we parameterize the function Vq rather than the function p, it will be 
useful to express the quantity e(xo), defined in Eq. ( |4.9|) and afterwards entering 
the important Eq. (4.18), also in terms of Vq. Obviously, using Eqs. (4.12) and 
Q4.34D one gets 



e(x )= / p(e,x )ede = - / V£(0t(t 



M 

Y 



d£ (4.40) 
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M 



■' 



Now, we can our results shortly summarize. If there are given some values of 



the free parameters c, 



q, < 



, s, ai, a2, then the corresponding proton and neutron 



s truct ure functions can be directly calculated accordin. 
(4.32), where the distribution /iy is given by Eq. (4.3E 



to Eqs. (|4.17|) , flUf) , 

the function e(xo) by 
Eq. ( [4.40D with the use of Eqs. ( [4.34|) , (|4.35| ), ( |4.39| ) and the functions <p, ip 2 



are: 



tp{x,x ) = 2x / T](£,x )d£, 



(4.41) 



V>i(x,Xo) = x / ??(£,Xo) 



x£ 



e + 4-x^ 



(4.42) 



where 



ip 2 (x,x ) = — r){£,x )—= =dg, 
2 A +^0-^) 



»7(f> zo) = c norm 0(i - Xq) e2 [ ^ 



e+xt 



(4.43) 



(4.44) 



The last expression is calculated from Eqs. ( 4.34| ), (4.35) and (4.39) with the 
use of the step function 0(x) =+!(—!) for x > (x < 0). 



5 Comparison with the experimental data 

Now we shall try to compare our formulas for the structure functions with the 
existing data. We shall not attempt to make a consistent, global fit of the free 
parameters based on some rigorous fitting procedure, but only show the set of 
optimal parameters obtained by their tentative varying on the computer "by 
hand" . Moreover, our constraint will be only agreement with the proton struc- 
ture functions F 2 and g%. It means that the parameter a, controlling asymmetry 
between the proton and neutron spin functions, must be somehow fixed before 
the fitting. For the first approximation we use the ST/ (6) value, a = 2/3 [see 
item 6c) in the previous section] . 

For a comparison with F 2 we use the parameterizations of the world data 
suggested in Q and pi, both taken for Q 2 = \QGeV 2 /c 2 . The data for g x are 



taken over from the paper 25 of the SMC Collaboration. After some checking 



on the computer, the optimal set of the free parameters is considered: 

c q = 0.43, r = -0.49, s = 6.5, a x = 1.6, a 2 = 1.5 (5.1) 
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Figure 6: Proton spin structure function gi(x) at Q 2 — lOGeV 2 / c 2 . The points 
represent experimental data [E5|, the curve is the result of our calculation. 
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Figure 7: Proton structure function F 2 (x) at Q 2 = lOGeV 2 /c 2 . The dotted and 
dashed curves represent the fits of the experimental data suggested in |2^] and 
[E5| . The full curve is the result of our calculation. 
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Figure 8: Proton structure function xg 2 (x) at Q 2 = lOGeV 2 / c 2 . The points 
represent the term xg 2 vw (x) (see text), the curve is the result of our calculation. 



Results of the calculation of the proton structure functions g\ and F 2 with 
these parameters are shown in Figs. [| together with the data. Let us 
remark, the experimental points for g\ correspond to the values evolved in |25| to 
Q 2 = lOGeV 2 1 c 2 . In the error bars all the quoted errors (statistical, systematic 
and those due to uncertainty of QCD evolution) are combined. Obviously, the 
agreement with the experimental data in both the figures can be considered 
very good, particularly if we take into account that our parameterization of 
the unknown distributions is perhaps the simplest possible and moreover, the 
parameters ( |5.lD still may not be optimal. 

Now, having "tuned" the free parameters by the g\ and F 2 , one can predict 
the remaining functions g 2 and F 2va i- The results are shown in Figs. |^, |9|. Our 
xg 2 surely does not contradict the experimental data ]2r| , which are compatible 
with zero - with statistical errors bigger, than the vertical range of the figure. 
Thus instead of the data, the comparison is done with Wandzura Wilczek p7| 
twist-2 term for xg\ vw , which is evaluated in [ p6| from the corresponding g\. It 
is obvious, that two completely different approaches give at least qualitatively 
very similar results. The proton valence function F 2va i in Fig. ^ is compared 
with the corresponding combination of the distributions xuv{x) and xdy(x) 
obtained (for Q 2 = 4GeV 2 / c 2 ) in |^8| by the standard global analysis: 

F 2va i(x) = ^xu v {x) + -xdv(x), J u v {x)dx = J dy(x)dx = 1. (5.2) 

Apparently, the agreement can be considered good. One can note, that the two 
different procedures, the standard one (uses input on F 2 , F£ N + QCD) and ours 
(uses input on F 2 , g± + our model) give a very similar picture of the function 
F 2va i(x), which is not directly measurable. 
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Figure 9: Proton structure function F2 Va i{x). The dashed curve represents the 
function based on the standard global analysis according to the relations (5.2). 
The full curve is the result of our calculation. 



6 Discussion 



Let us make a few comments on the obtained results. First of all, it should 
be pointed out, that our structure functions in Figs. are calculated on the 
basis of very simple parameterization of the unknown distributions /j,(xq) and 
Vq(x,Xo), but on the other hand it is essential, that the contributions from the 
individual components of the quark distribution correctly take into account the 
intrinsic quark motion, which is particularly important for the spin structure 
function. The effect of this motion on g\ is demonstrated in Fig. |^ and the 
fact, that we succeeded to achieve a good agreement with the data in Fig. |^ 
is just thanks to this effect. For a better insight, how our structure functions 
are generated, in the Fig. ^ we have displayed the initial distribution function 
Vq{x,Xq) drawn for a few values Xo, together with the corresponding structure 
functions F2, <7i, xg^- The complete structure functions are their superpositions 
- weighted by the corresponding way with the use of the distribution /j,v{xq). 

Further, also some other assumptions of the model are possibly oversim- 
plified, for a more precise calculation, at least some of them could be rightly 
modified - but at a price of introducing the additional free parameters. For 
example, the constant w sea [see Eq. (4.14)] should take into account some sup- 
pression of the s— quarks p8[ and probably should allow a weak dependence 
on xq. Also for the constant c q [see Eq. (4.6)] some xq~ dependence should 
be allowed. Concerning this constant, let us make one more comment. The 
standard global fit [^ suggests (at Q 2 = lOGeV 2 /c 2 ) the quarks carry ~ 56% 
of the nucleon energy and our fitted value c q from the Eq. (|4.6| ) is roughly 43%. 
This difference is mainly due to the different relations between the distribution 
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Figure 10: Distribution functions Vq(x,xq) drawn for x a = 
0.005,0.015,0.05,0.15,0.5 and the structure functions generated corre- 
spondingly. The calculation is based on the Eqs. ( [4.34 ), (4.12), ( 4.30| ) and 
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and structure function in both the approaches, see Eqs. (3.38), ( [3.62 ). The 



second relation (valid for a subset of quarks with effective mass xq), multiplied 
by x 2 and then integrated by parts gives 

1 1 r 1 fx 2 
1 / ip i \ { o , x o 



xF(x,xo)dx = - F 2 (x,xo) 3H — ~ ) dx, (6.1) 

4 Jxl \ xI J 

which for the static quarks [F(x,Xq) ~ F 2 (x,Xq)/x ~ 6(x — a; ), see discussion 
after Eq. ( 3.62| )1 coincides with the standard relation. Nevertheless, generally 



both the relations imply different rate of the nucleon energy carried by quarks. 
One can check numerically that for our F 2 (x,xq) in a dominant region of xo 
the term (x/xo) 2 in the integral (ELI) plays a minor role (see also Fig. [To| , 



positions of the maxima of F 2 's are above the corresponding xo, in particular 
for lower xq), so as a result we get 3/4 of the standard estimation of the quark 
contribution to the nucleon energy. This ratio agrees with the ratio obtained 
from the corresponding fits: 3/4 ~ 43% / 56%. 

In the previous section we have assumed the same shape for the valence terms 
related to the u and d— quarks. This assumption together with our premise a = 



2/3 (i.e. SU(6) symmetry) in an accordance with Eq. (4.24) give H n (po) = 
and g™ (x) = correspondingly. On the other hand it is known, that the neutron 
structure function <7™(cc) ^ 0, even if |r™| is substantially less than T^. A more 
consistent approach (a = 2/3 but d ^ u) would give 

H n ( Po ) = ^(-u n ( Po ) + d"(po)) = ^(-d"(po) + « P M), (6.2) 

therefore d ^ u implies <7™ (x) ^ 0. Actually, the global fit analysis proves that 
d^ al (x) is slightly "narrower" than u p val {x). It means, considering qualitatively, 
in accordance with the equation above in the function g™ (x) the negative term 
should dominate for smaller x, which does not contradict the data. A proper 
accounting for this difference into the model should enable to calculate con- 
sistently in a better approximation not only the proton and neutron structure 
functions F 2 ,gi,g2, but also the neutrino structure functions. Apparently then 
one could make a "super-global" fit covering the both unpolarized and polar- 
ized DIS data. As a result, the flavor-dependent quark distributions Vo(x,Xq) 
[or equivalcntly p(e, xq)] together with the corresponding effective mass distri- 
butions and the parameter a controlling the relative spin contribution of the u— 
and d— quarks, could be obtained. 

Finally, let us point out, inclusion the spin structure function into the fit 
in our model enables to obtain some information about the distribution of the 
quark effective masses. Within our approach there are two distributions, /zy 
and (i, relevant for the description of the quark effective masses in the nucleon. 
The extrapolation of our parameterization for the fi distribution with the use 
of the relations a w mj (E kln ) and (fd||) , jO^) , (^39|) , fTi"!) give for x -> : 



K x o) _^ x o y J£q 

e(x ) Mx + (E kin ) ^ jlnxol 1 



M*o) ~ ^ - „_ ; - ^-ct, (6.3) 
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which implies the extrapolated /i is not integrable in this limit. On the other 



hand, the basic distribution /iy, parameterized by Eq. (4.33) with the r, s from 
the set flj.ip and with the use of the known relation zT(z) = T(z + 1) can give 
an estimate of the mean value: 

^°)v = 1 9 - 0-064, (6.4) 

i.e. (m) ~ QQMeV for Q 2 = 10GeV 2 /c 2 . The corresponding kinetic term 
calculated as 



(E k in) v = J nv(x )(e(x ) - Mx )dx (6.5) 

gives a similar number (~ 60-MeV). This number agree very well with the 
corresponding temperature obtained in the statistical model (2^|. The In- 
dependence is involved only in the distribution /iy (xq, Q 2 ), i.e. in our param- 



eterization (4.33) only via the powers r(Q 2 ), s(Q 2 ). It follows, the structure 
functions, which enhance in a low— a; region for increasing Q 2 , must be gen- 
erated by the distribution hv{xq,Q 2 ) in which the mean effective mass (xq) v 
drops for increasing Q 2 - in the qualitative agreement with an intuitive expec- 
tation. 



7 Summary 

In the present paper we proposed an alternative, covariant formulation of the 
QPM. The initial postulates of the standard and our approach are basically 
the same, despite of that the relations between the structure and distribution 
functions obtained in both the approaches are not identical. It is due to the 
fact, that in the standard approach the intrinsic quark motion is effectively 
suppressed by the use of the approximation = xP^. On the other hand, 
we have shown the master equations can be solved without the use of the this 
approximation, so in the corresponding solution the quark intrinsic motion is 
consistently taken into account. In particular, we have suggested, that the quark 
intrinsic motion can substantially reduce the structure function g\. 

On the basis of the obtained relations (a priori valid for the version of naive 
covariant QPM - with nonstatic quarks on mass shell) we propose the model, in 
which the distributions (/zy,Vb) reflecting the parton dynamics are introduced 
with some free parameters. With the use of this model we calculated simultane- 
ously the proton structure functions i*2, F 2va i, 9i, gi-, assuming only the valence 
quarks term contributes to the proton spin. Then by a comparison with the 
data (F 2 ,g 1 ;Q 2 = lOGeV 2 / c 2 ) we fixed the free parameters. We found out: 

1) Both the unpolarized structure functions can be well reproduced by the 
model. The comparison is done with the data on F 2 and with the F 2va i obtained 
from the standard global analysis data. 

2) At the same time, the model well agrees with the data on gi and the 
calculated g 2 does not contradict the existing experimental data. 
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3) Analysis of the fixed parameters within our approach suggests: 

i) The quarks carry less the proton energy (almost by the factor 3/4), than 
estimated from the standard analysis. 

ii) The average effective mass related to the valence quark term can be roughly 
60MeV and a similar energy can be ascribed to the corresponding motion. 
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